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1 INTRODUCTION 



The theoretical question which will be addressed in this talk is the computation 
of quark and gluon distribution functions at small Bjorken x. The reason why 
such computations may be possible is because at small values of Bjorken x, which 
correspond to central values of rapidity, the density of partons per unit area satisfies 

This density of partons per unit area is the only local parameter with dimensions 
of an energy squared which we can construct, and the coupling therefore should be 
evaluated at this scale. Defining 

.9 1 diV 

the strong coupling constant at this scale must be 

a,(A) « 1 (3) 

so that it may be possible to formulate the problems in weak coupling. 

There are several problems which might be solved if one can compute the quark 
and gluon distribution functions at small Bjorken x. At small x, it is expected that 
the gluon distribution function for a single proton behaves as 

~ (4) 

This implies that the distribution of partons for a proton diverges as x — > 0. If this 
is the case, then this problem should be possible to analyze using weak coupling 
methods. The distribution functions at small x have recently been measured at 
HERA, and do seem to be singular at small x. 

In addition to understanding the Lipatov enhancement, it would also be useful 
to compute quantities such as the ratio of sea quark to gluon distributions. For 
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example, if the typical energy scale has been increased due to a high density of 
partons, we should expect that heavier quarks will become of increasing importance. 
Charm quark production might for example become substantial. 

Another interesting physical problem is deep inelastic scattering and di-lepton 
production using nuclei with A» 1. If the value of Bjorken x « A~^/^, then the 
nucleus as seen by a parton moving with that value of x is Lorentz contracted to a 
scale size which is much smaller than the wavelength of the parton in a frame comov- 
ing with its longitudinal momentum. It is expected that in this kinematic region, 
there should be non-trivial effects which might screen the effects of the valence 
nuclear matter distribution. On the other hand, we expect that the distribution 
functions for partons should become large for large nuclei, and if the effects due to 
screening can be ignored (as we will see they can be when we do the computation) 
then 

(5) 

as A — > oo. In addition to the questions about the Lipatov enhancement and 
the ratio of sea quarks to glue, one can also ask about the A dependence of the 
distribution functions. 

Finally, there is the problem of determining the initial conditions for quarks and 
gluons in heavy ion collisions. If one can understand the distribution functions, then 
these may provide information about the boundary conditions for the evolution of 
the matter into a quark-gluon plasma. Recall that since the density of partons per 
unit area is the only scale in the problem, and it goes like A^^^, the energy density 
will have to scale as E/V ^ A"^/^, where this is the energy density scale at a time 
which corresponds to the dimensional scale constructed from the density of partons 
per unit area. In order for these considerations to be valid, the Lorentz contracted 
size of the nucleus must be smaller than the size scale constructed from the density 
of partons per unit area. This is Ecm/A » s/A which for large nuclei requires 
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that Ecm/A >> 50 GeV which is within the range accessible at RHIC. 

In what fohows, I will describe how to compute the quark and gluon distribution 
functions for very large nuclei. We have not yet succeeded in being able to compute 
for a single hadron. It may also be true that our weak coupling analysis may be 
marginal at best for realistic values of A, and therefore our results may only be 
useful as a theoretical laboratory which will give us some insight into the structure 
which we might expect for realistic nuclei. If it is true that weak coupling techniques 
cannot be used for the earliest stages of a heavy ion collision, then it also will imply 
that such techniques are probably never applicable in the subsequent evolution of 
the matter produced in such a collision. 



2 Summary of Results 

Before going into a discussion of how to analyze the problem of computing distri- 
bution functions for very large nuclei, I will first summarize our results to date.[§] 
First I define 

A j\jvalence 

f^' = - I-IA^/' Fm-' (6) 

3 vr/t"^ 

which is the the density of valence quark color charge squared per unit area. I will 
show that there is a many body theory which describes the quark and gluon distri- 
bution functions so long as we restrict our attention to parton transverse momenta 
which satisfy 

4 « /"^ (7) 
and we require that we are at small values of x 

X « (8) 

There are two expansion parameters for this many body theory as(/^) and as{n)lJ./qT 
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In the range as(//)/x << qt « n the gluon distribution function may be 
evaluated to lowest order in as{^Ji) as 

1 dN _ «,(iV,^-l) /x^ 
7ri?^ dxcPgr tt^ xqj^ 

This is just the Weizsacker- Williams distribution function for a Lorentz boosted 
distribution of Coulombic gluons, scaled by the factor /x^ which is the average value 
of the charge squared per unit area. 

The result above is however only a formal result if it is true that there is a 
Lipatov enhancement. The Lipatov enhancement is of the form x"*^"" If we formally 
expand this as a series in ag, we find the expansion parameter is as{iJ,)ln{l/x). 
Since we have assumed that x « A~^/^ and fj? ~ AV3, this expansion is not 
well behaved. To correct for this behavior if it occurs, it will be necessary to go 
beyond the naive weak coupling expansion. It should however still be true that 
weak coupling methods are still applicable. 

If we sum to all orders in but to first order in a^, we have shown that the 
gluon distribution function is of the form 

- —H{q^/ain') (10) 



■kE? dxd'^QT TT^ X as 

where the function liniy-^ooH {y) — ^ 1/y The function is a correlation function for 
an ultraviolet finite two dimensional Euclidean quantum field theory. The strong 
coupling limit of this field theory is the small qr limit. In this limit, we expect 
disorder since the theory is strongly coupled. This implies exponentially falling 
correlations in coordinate space which implies that H{Q) should be finite. 

3 Setting the Problem Up 

To begin computing the distribution functions, we will assume we are working in 
a frame where the nucleus is moving close to the speed of light. In this frame the 
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nucleus appears as a Lorentz contracted pancake. The value of Bjorkcn x in this 
frame is essentially the ratio of longitudinal momentum of the parton to that of the 
projectile nucleus per nucleon. 

Since the variation in the nuclear valence quark distribution function with re- 
spect to transverse coordinate is 

1 dN 1 



(11) 



N dr Rnuc 

and since the average scale of transverse momenta is qt » 1/Rnuc-, we see that 
locally, the variation in the transverse nuclear matter distribution can be ignored. 
Therefore to compute the local properties of the gluon distributions function as 
a function of transverse coordinate,we need only consider a problem where the 
transverse matter distribution is assumed to be uniform and of infinite extent in 
the transverse direction. The problem is therefore of a sheet of valence quarks 
uniformly distributed on a thin sheet of infinite extent in the transverse direction. 
We will take the density of valence quarks to be N quark /t^R^ ~ A^/^ . 

The natural variable with which to analyze the dynamics are light cone vari- 
ables, 

° = <'^) 

If we let A;+ be a parton light cone momentum and p+ that of a projectile nu- 
cleon, then X = Ip^ In these variables, instead of constructing eigenstates of 
the Hamiltonian, it is simpler to construct the eigenstates of the generator of x"*" 
transformations , 

P- = ±(H-P^) (13) 

Finally, when constructing the light cone Hamiltonian and action for the theory 
in terms of these variables, it is simplest to work in light cone gauge, 

= -A+ = (14) 
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In this gauge, the hght cone Hamiltonian has the form 

P- = j Sx^F^ + ^{pf + Dt-Et)^{pf + Dt-Et) 



In this equation, 



+ ^V^(M-7-Pt);^(M + 7-Pt)V' (15) 



PM = i^M -gr-A*" (16) 
i 



and Ek = —d-Ak- The Lagrangean is generated in the usual way by adding in 
—itjj^d+il! + Ekd+Ak The fermion charge density is 

= V^^+r'^V (17) 

The quantized fields are 

Jk+>o {^T^r 



r H^k X 

A^x) = (<(fc)e-^- + af(fc)e-^-) (18) 

Jk+>o {2-!ryV2k+ 

The last ingredient we need to set the problem up is how to describe the valence 
quarks. These quarks are travehng close to the speed of light. For small values of 
as, these quarks occasionally emit a small x gluon. This emission does not change 

the path of the valence quark. It should therefore be a good approximation to treat 
the trajectories for the valence quarks as straight line propagation at light velocity, 
that is 

J+ = Pa{x+,XT)S{x-) (19) 

Unlike the case for QED, in general in QCD the charge density will have to depend 
on the time a;"*" since the extended current conservation condition requires that 

+ r^^AbQc = (20) 

This forces the charge to rotate as 

T-Q{x+) = U{x+)t ■ Q{0)U^ {x+) (21) 



where 

U{x+) = Texp dy+ r • A+{y+)^ (22) 

To summarize, we have shown the problem which must be solved is to compute 
the ground state expectation values for a system with the valence quarks traveling 
with light velocity localized along an infinite sheet in the transverse space. The 
density of the valence quarks is uniform. This problem is well posed since the con- 
straint on the valence quarks is equivalent to specifying the space-time coordinates 
of the electromagnetic charge and baryon number. These operators commute with 
the QCD Hamiltonian. 

4 The Example of QED 

In QED, the problem we want to solve is the photon structure function generated 
by a fast moving electron. We will ignore pair production of electron-positron pairs. 
The source for the electron is x"*" independent and is 

Pe = e6ix-)S^^\xT) (23) 

The light cone Hamiltonian is 

P- = J d?x QfI + i(p, + Vt • ET)^{pe + Vt • Et)) (24) 

For this Hamiltonian, the ground state is a coherent state 

I * >= Cexp J (fxA°P{x)E''\x)j \ > (25) 

Letting operate on this state, we see that the ground state has P~ = so that 
the classical field is purely longitudinal and 

Bt = 

Vt-Et = -epe (26) 
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The solution for these equations are that 

In space-time, the vector potential is ^(a;~)VA which is for a;~ < vanishing and a 
pure gauge for x~ > 

The field above is precisely the Weizsacker- Williams field for the Lorentz boosted 
Coulomb field. The distribution function for the photons can be computed and is 

1 

(2V)3 



F{k+,kT) = <a\k)a{k)> 



2e2 1 

(28) 



(27r)3 k+k^ 



or 

a 1 



^(^'^*) = ^^ (29) 

5 The Distribution Functions for QCD 

For QCD, we have not been successful in constructing the ground state wavefunction 
in the presence of the external source corresponding to the valence quarks. We have 
however concentrated on computing ground state expectation values. To do this, 
we first note that 

Z = lirriT^oo ^<N\ e*^^" | N > (30) 

N 

will project onto the ground state. The sum over N here includes a sum over the 
color labels of the external source of color charge generated by the valence quarks. 

To do the sum over N, wc break our transverse space into a grid of squares 
with size cPx >> ttR"^ /Ng^^^ark ~ A~^/^Fm'^ Our approximation will therefore only 
be good when we look at transverse momentum scales where « In these 
limits, the number of valence quarks in each square is much larger than 1. If this is 
the case, then typically the charge in each square will be much larger than 1. If Q 
is this charge, then >> Q ~ [Q, Q] so that the charge may be treated classically. 



If the total charge of interest is also much less than the maximum possible 
charge in the square, then the density of states for charge Q is e~*^^/^'^^. Summing 
over the states in the definition of Z is therefore equivalent to inserting into the 
path integral the integration 

where p is the charge density per unit area. 

For such a Gaussian charge distribution we have 

< P%x)p\y) >= p^S-'S^^^ {XT - yr) (32) 

where it is straightforward to estimate p? = 1.1 A^^^ Fm? 

The problem we must solve is therefore that described by the theory in the 
presence of an arbitrary external source of surface charge on the light cone and then 
integrating over all possible values of the charge. The fields generated correspond 
to a stochastic source of charge. The problem has therefore been reduced to a many 
body theory. It is possible to integrate out the sources entirely and get an action 
in terms of the quark and gluon fields with a term involving /x^ which is associated 
with the valence quark charge density. The value of fi^ sets the scale for the coupling 
constant. For large fi^ the coupling is small and weak coupling methods should be 
reliable 

We can now compute the gluon distribution function to lowest order in and 
to lowest order in p?. We first compute the change in the propagator induced by 
the sources. This is 

S < A^A^ >== J [dp] < P"(x) p\x) > (33) 



or 
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Using this form of the propagator, it is now easy to show that 

This is just the Weizsacker-Wilhams distribution function weighted by the average 
charge squared per unit area. This result reflect the RMS fluctuations in the stochas- 
tic background field induced by the source associated with the valence quarks. 

We will soon see that this result is valid only in the range of momentum where 
alii'^ « « ji"^ The last term in the previous limit is just the region of validity 
of our derivation of the many body theory. The first term is the limit of validity of 
assuming that fi^ is small and that one can expand to first order in this quantity. 

We can generalize our results to all orders in a^/x and first order in as To do 
this, we solve the classical problem of computing the fields in the background of an 
arbitrary source and then integrate over the source. This classical problem will be 
accurate to first order in ag 

We must solve the equations of motion 

D^F''^ = g5^+5{x-)p{x+,XT) (36) 

There is a solution of these equations of motions with = A- = which also has 
Ft = 0. This solution is of the form 

Aj{x-^XT) = e{x~)aj (37) 

The condition that Ft = is equivalent to the condition that the field aj is a gauge 
transformation of the vacuum configuration of a two dimensional gauge theory. The 
condition that Dt ■ Et = —gp is equivalent to the two dimensional gauge condition 

V-a = -gp (38) 

The above field configuration may be written as 

T-a = -—U(xTWiUUxT) (39) 
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and the gauge condition is 

V{-UVU^) = ig^p (40) 
The integration over the sources may be written as 

j[dU]exp(^--^tr(vT-U-^VTU'^y^ (41) 

There is of course a Fadeev-Popov determinant for this measure which comes from 
restricting to Feynman gauge, but we wiU not be concerned with this measure 
here. (The determinant does not affect the arguments we present for the vaHdity of 
perturbative expansions nor the scahng behavior in kx-) 

Note that the couphng constant for this theory is g^jj, so that the expansion pa- 
rameter is g^ji/qx- The fluctuations over the external charge generate an ultraviolet 
finite theory. The correlation function of UVU^ /g generates the modifications due 
to the sources of the gluon propagator. This correlation function should die expo- 
nentially at long distances corresponding to small coupling , that is, in momentum 
space, the correlation function should be finite at zero momentum. 

The step function 0{x~) in the solution of the classical field equations guarantees 
that the distribution of gluons is proportional to when the theory is solved to 
first order in ag The gluon distribution function is therefore of the form 

1 dN iV2-l 1 1 



ttB? dxdPkx TT^ X as 
where limy^ooH{y) = 1/y 



- -H{ki/ay) (42) 



6 Summary 

There are of course many problems which must be solved before this approach 
can provide a realistic theory of the distribution functions at small x. The first 
order corrections in must at least be understood. This will generate the induced 
contribution of sea quarks. It will also presumably lead to the Lipatov enhancement. 
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This must be nontrivial since for this theory where the valence quarks are locahzed 
to a delta function on the light cone, there is no scale of P"*", the momentum per 
nucleon of the nucleon in the nucleus. This can only arise from the cutoff dependence 
of a regularized delta function. Such a dependance will not affect the physics to 
leading order in as nevertheless, since if we change the cutoff by a finite amount 
(A/A;+)*^"^^ changes by only an amount proportional to ag 

In addition to the above problems, the issue of how large a transverse mo- 
mentum one can use to define a sensible theory remains open. Even though our 
approximations were only valid for qr « ^ the validity of assuming straight line 
trajectories along the light cone should be valid in a much broader region. The ex- 
tent to which the theory can be solved in the large momentum transfer region and 
the extent of the region of validity of the external source approximation remains 
open. 

Finally, there is the issue of actually computing structure functions for deep 
inelastic scattering or Drell-Yan particle production. We have here only discussed 
computing the expectation values of quark distributions in the ground state wave- 
function. When a probe of momentum is introduced, in general there will be two 
parameters of interest and /x^. For q^ » ji^, presumably one can analyze the 
distribution functions using Altarelli-Parisi equations. For q^ < the situation 
must be more complicated since here the momentum of the quarks inside the hadron 
wavefunction are important and one is not in the scaling region. 
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